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Q ! Abstract 



We have investigated the thermodynamical properties of the Universe with dark energy. Adopt- 
ing the usual assumption in deriving the constant co-moving entropy density that the physical 



volume and the temperature are independent, we observed some strange thermodynamical behav- 

t-h ; 

\ iors. However, these strange behaviors disappeared if we consider the realistic situation that the 

physical volume and the temperature of the Universe are related. Based on the well known corre- 
spondence between the Friedmann equation and the first law of thermodynamics of the apparent 
horizon, we argued that the apparent horizon is the physical horizon in dealing with thermodynam- 
ics problems. We have concentrated on the volume of the Universe within the apparent horizon 
and considered that the Universe is in thermal equilibrium with the Hawking temperature on the 
apparent horizon. For dark energy with w > — 1, the holographic principle and the generalized 
second law are always respected. 
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I. INTRODUCTION 



Since the great discovery that the Universe is experiencing accelerated expansion driven 
by dark energy (DE) from the type la supernova (SN la) observations in 1998 1], many 
works have been done in pursuit of understanding this spectacular phenomena. The later 
more accurate SN la data [2|, together with the Wilkinson Microwave Anisotropy Probe 
data [3] and the Sloan Digital Sky Survey data [4J indicate that the Universe is almost 
spatially flat and further support the existence of DE which contributes about 72% of the 
matter content of the present universe within the framework of Einstein's general relativ- 
ity. Despite the robust observational evidence for the existence of DE, DE is still a major 
puzzle of modern cosmology. Except knowing that DE has negative pressure, we know little 
about its theoretical nature and origin. In this work we would like to study the DE from 
thermodynamical considerations. 

In the semiclassical quantum description of black hole physics, it was found that black 

holes emit Hawking radiation with a temperature proportional to their surface gravity at 

the event horizon and they have an entropy which is one quarter of the area of the event 

horizon in Planck unit [5L The temperature, entropy and mass of black holes satisfy the first 

I 

law of thermodynamics [6( . On the other hand, it was shown that the Einstein equation can 
be derived from the first law of thermodynamics by assuming the proportionality of entropy 
and the horizon The Einstein equation for the nonlinear gravitational theory 

f{R) was also derived from the first law of thermodynamics with some non-equilibrium 
corrections [8] . For a general static spherically symmetric space-time, Padmanabhan showed 
that the Einstein equation at the horizon gives the first law of thermodynamics on the 
horizon [9j. The study on the relation between the Einstein equation and the first law 
of thermodynamics has been generalized to the cosmological context where it was shown 
that the first law of thermodynamics on the apparent horizon ta can be derived from the 
Friedmann equation and vice versa if we take the Hawking temperature Ta = l/27rf J 4 and the 
entropy Sa = 7rr^/ G on the apparent horizon [lO]. Furthermore, the equivalence between the 
first law of thermodynamics and Friedmann equation was also found for gravity with Gauss- 



Bonnet term and the Lovelock gravity theory [10|, llll]. These results disclosed that there 
may be a deep relationship between the thermodynamics of the horizon and the Einstein 
equation, which might shed some light on the properties of DE. 
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Motivated by the black hole thermodynamics, Bekenstein proposed a universal entropy 
bound for a weakly self-gravitating physical system in 1981 [12j. 't Hooft and Susskind 
subsequently built an influential holographic principle 13J relating the maximum number 



of degrees of freedom in a volume to its boundary surface area 



14 



15) . The extension of 



the holographicprinciple to a general cosmological setting was first addressed by Fischler 



and Susskind 



161 ] and subsequently got modified by many authors [17], [l8|. The idea of 



holography is viewed as a real conceptual change in our thinking about gravity [19]. There 
have been some examples of applying the holographic principle to understand cosmological 



problems. It is interesting to note that ho 
ical constant in a large class of universes 



ography implies a possible value of the cosmolog- 



201 ] . In an inhomogeneous cosmology holography 



2l|. The idea of 



was also realized as a useful tool to select physically acceptable models 
holography has further been applied to the study of inflation and gives possible upper limits 
to the number of e- folds 22J. Recently, holography has again been employed to investigate 



the DE 23J and it is expected that the holographic principle could serve as an effective way 
to help us understand the DE. 

The Universe can be considered as a thermodynamical system. The thermodynamics in 



de Sitter space-time was first investigated by Gibbons and Hawking in [241 ] . In a spatially flat 
de Sitter space-time, the event horizon and the apparent horizon of the Universe coincide and 
there is only one cosmological horizon. It was found that the area of the cosmological horizon 
can be interpreted as the entropy or lack of information of the observer about the regions 
which he or she cannot see, and an observer with a particle detector will observe a background 
of thermal radiation coming from the cosmological horizon. The thermodynamical study of 
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27| . When the 



the Universe has been extended to the quasi-de Sitter space in [25, 
apparent horizon and the event horizon of the Universe are different, it was found that the 
first law and the second law of thermodynamics hold on the apparent horizon, while break 



down if we consider the event horizon 



271 ] . Due to this subtlety, we will concentrate our 



attention on the thermodynamical properties of the Universe enveloped by the apparent 
horizon. We will derive the temperature and entropy of the matter contents inside the 
apparent horizon from the first law of thermodynamics and discuss the holographic entropy 
bound and the generalized second law (GSL) of thermodynamics for the Universe with DE. 

In the discussion of the thermodynamical properties of the matter contents, it is usually 
assumed that the physical volume and temperature of the Universe are independent and 
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by using the integrability condition d 2 S/dTdV = d 2 S/dVdT one has the constant co- 
moving entropy density a and the relation Ts = p + p for the physical entropy density s. It 
follows that either the entropy or the temperature is negative for the phantom. Applying the 
constant co-moving entropy density to the first law of thermodynamics for DE with constant 
equation of state w = p/p, the thermodynamical properties of DE and the thermal spectral 



distribution were discussed in 



28 



. The thermodynamics of the DE with constant w in the 



range — 1 < w — 1/3 was also discussed in 29j, and the thermodynamics of phantom with 
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321 • 



w < — 1 was investigated in [30]. Some other discussions on this topic can be found in 
However, in general, when we consider the thermal equilibrium state of the Universe, the 
temperature of the Universe is associated with the horizon, then the integrability condition 
cannot be used to derive the constancy of the co-moving entropy density from the first law 
of thermodynamics. Some assumptions on the temperature or entropy are needed to derive 
the thermodynamical properties of DE. With the requirement that the entropy of DE is 
bounded from above and the assumption that it could not increase faster than that of cosmic 
microwave background radiation, we will investigate the thermodynamical properties of DE 
in the thermal equilibrium universe. Employing the holographic entropy bound, we will 
also discuss the GSL for the Universe with DE. The study of GSL in a phantom dominated 
universe was carried out in [33j . We will address the thermodynamics of DE by considering 
the DE models with constant w and the generalized Chaplygin gas (GCG) model which 
extrapolates between a dust and a cosmo logical constant 34]. 

The paper is organized as follows. In section [Til we review the equivalence between 
Friedmann equation and the first law of thermodynamics of the apparent horizon. In section 
IHIt we discuss the thermodynamical consequence by assuming the physical volume and 
temperature are independent. In section HVl considering the physical volume is a function of 
the temperature, we investigate the thermodynamical properties for the DE with constant w 
and the generalized Chaplygin gas. We discuss the GSL of thermodynamics for the Universe 
with DE in section |VJ The summary and discussion will be presented in the last section IVI1 

II. FIRST LAW OF THERMODYNAMICS ON THE APPARENT HORIZON 

For a spherically symmetric space-time with the metric ds 2 = g a bdx a dx b + r 2 dQ 2 , where 
the unit spherical metric dVL 2 = d6 2 + sin 2 ^^ 2 , we define the Misner-Sharp mass M as 



4 



351,136] 

j*W = f = l-^. (1) 

r 

The mass M defined in Eq. (OQ) was shown to be the active gravitational ener gy in the 
vacuum, small sphere, large sphere, Newtonian, test particle and relativistic limits 37j]. The 
Einstein field equations tell us [35] 

M ;a = 47rf 2 (T a 6 - 8 b X)r., b , (2) 

where the semicolon denotes covariant derivative with respect to the two dimensional metric 
g a t). The apparent horizon is determined from / = 0. We define the dynamic surface gravity 
at the apparent horizon as 

(3) 



So the Hawking temperature at the apparent horizon is T = \k,\/2tt. The surface gravity 
k defined in Eq. ([3]) has the same form as the Newtonian surface gravity. In static black 
hole physics, the horizon is located at f(f) = 1 — 2GM/r = which gives f = 2GM. 
The surface gravity is found to be k = —f'/2 = — 1/AGM and the Hawking temperature 
is T = \/87rGM. For the de Sitter universe, f(r) = 1 — H 2 r 2 , the event horizon which 
coincides with the apparent horizon is ta = H^ 1 . The surface gravity at the apparent 
(event) horizon is found to be k — —f'/2 = f^ 1 and the Hawking temperature T = l/2nfA- 
For static spherically symmetric space-time, it was shown that Einstein equation at the 
apparent horizon fj gives the first law of thermodynamics of the apparent horizon with 
the identification of the temperature as T = |/ / (f i )|/47r and the entropy as S = A/ AG, 
where A = Anrf is the area of the horizon [9]. It was shown that the apparent horizon is a 
good boundary holding thermodynamical laws in the Universe driven by quintessence and 
phantom 271 ]. 

For the Friedmann- Robert son- Walker (FRW) metric, we have 




9ab = I _ n _ . nv I , r = ar. (4) 

So 



The apparent horizon is 



f = g ab r, a r, = l-[H 2 + \)r 2 . (5) 



The surface gravity at the apparent horizon is k = —f/2 = l/r A and the Hawking temper- 
ature is 

Ta = (7) 
2nr A 

By using the FRW metric (jlj), the definition of the Misner-Sharp mass (JTJ), and the perfect 
fluid plus the cosmological constant for the matter, we get the Friedmann equations from 
the mass formulae (J2J), 

H* + - = — P + A, (8) 
p + 3H(p + p) = 0, (9) 

H- \ = -47rG(p + p). (10) 

or 

Eq. (TTOl) can be derived from Eqs. (ED and ([9]), so there are only two independent equations. 
Combining Eqs. (EI) and ( jTOi) . we have 

f A = -r\H(H - k/a 2 ) = AnGr\H{p + p). (11) 

The energy flow through the apparent horizon is 

-dE= An? A T^k»k v dt = kirr\{p + p)Hdt, (12) 

where the (approximate) Killing vector or the (approximate) generators of the horizon, the 
future directed ingoing null vector field k^ = (1, —Hr, 0, 0). Note that dE = k a M. a with 
k a r. a = 0. Therefore, we get the first law of thermodynamics of the apparent horizon 

1 A / Airr 2 \ 

T A dS A = — - — ^ = f A /Gdt = -dE, (13) 



2%f A dt V AG 

if we define the entropy of the apparent horizon as S A = A/ AG = 7rf A /G. It is clear that the 
first law of thermodynamics can be derived from Einstein equations. More interestingly, Ein- 
stein equations can be obtained from the first law. By using the first law of thermodynamics, 
we get 

r A = 4ttG? a H(p + p) = -f\H{H - k/a 2 ). (14) 

So 

H — k/a 2 = —AnG(p + p). (15) 

Combining Eq. (|T5|) with the energy conservation equation (Q, we can derive Friedmann 
equation (jSJ). In this derivation, there is also an integration constant which is set to be the 



cosmological constant. The derivation of Friedmann equation from thermodynamics was 



also shown for gravity with Gauss-Bonnet term and Lovelock gravity theories The 
equivalence between the Friedmann equation and the first law of thermodynamics of the 
apparent horizon for other gravitational theories, like the higher order gravity f(R) and the 
Brans-Dicke theory, remains to be an open problem. 



There are other general discussions on the first law of thermodynamics. In 37], the 
author derived a so called unified first law from the mass formulae (JTJ), 

V a M = A^ a + WV a V, (16) 

where the energy density W = — ~T", the energy-supply vector ty a = T^dbf + W d a f : , and the 
areal volume is determined by W = AVr. The first term in the unified first law may 
be interpreted as an energy-supply term, analogous to the heat-supply term in the classical 
first law of thermodynamics, while the second term VTVU may be interpreted as a work 
term. If we define the variation along the trapping horizon as 5M = z a, VaM, where z is a 
vector tangent to the trapping horizon and the normalization of z is irrelevant, then it was 
shown that Az a ^ a = Ki5A/8ir and the unified first law becomes 

5M = — 6A + W6V, (17) 
8% 

where K\ = Or/2 = M/r 2 — AnrW. For cosmological models with FRW metric, K\ ^ k 
and the above first law (TT7T) is different from the one (fTBl derived from the Friedmann 
equation directly. Although the unified first law seems to be true for more general geometry 
with trapping horizon, our definitions of the surface gravity and the Hawking temperature 
(J3]) are appropriate for use in the cosmological context since the first law built reflects the 
requirement of the Einstein field equation. 

Since at the apparent horizon, / = 0, so the vector V a / is normal to the apparent horizon, 
while the orthogonal vector ( a = e ab Vbf is tangent to it. If we define the variation at the 



apparent horizon as 5M = £ a V a M, then we have 



26) 



5E = ( a M, a = %~r A C~r A , a = ^—C{7tf\), a = ^6S. (18) 
r A 4-kGta 2 



In the above derivation, only the definitions of M and f A are used. Although the formulae 
( I18p looks in form like the first law of thermodynamics, it is different from Eq. ( TTBl . The 
physical meanings of 5E and 5S are different from dE and dS before. Further, Eq. (fl8|) 
cannot reflect the dynamical property of the Universe. 



Thus, in the cosmological context, the apparent horizon is the physical proper size in 
discussing thermodynamics and the appropriate description of the first law on the apparent 
horizon is expressed in Eq. (I13p with the definitions of the temperature in Eq. (j7j) and the 
entropy S = nr\/G on the apparent horizon. The first law for more general spherically 
symmetric space-times in other gravity theories needs further investigation. 



III. CONSTANT CO-MOVING ENTROPY DENSITY 

In an expanding universe, the first law of thermodynamics is 

TdS(T, V) = Td(sV) = d(p(T)V) + p(T)dV, (19) 

where s is the entropy density and V is the physical volume. The energy density and 
the pressure are assumed to be functions of temperature only, and the volume V and 
the temperature T are usually assumed to be independent. The integrability condition 
d 2 S/dVdT = d 2 S/dTdV relates the energy density and pressure 

dT~ T ■ W 

For the cosmological constant p = —p, the first law implies dS = d[V(p + p)/T] = 0. 
Substitute Eq. ([20]) into Eq. ([19]), we get 

a = P -±?, for p ± -p. (21) 

For phantom field, we see that either the temperature or the entropy is negative. Combining 
Eqs. (im ([20]) and (l2Tfl . we get 

dp dp 

ds = — = s^—, for p ^ -p. (22) 
T p + p 

Therefore, once we are given the equation of state p/p, we can get the relationship between 
the entropy density s and the energy density p, thence we can find the relationship between 
the temperature and the energy density. If we substitute the energy conservation equation 
([9]) into Eq. ([22]), we get 

cr = sa 3 = constant. (23) 

Therefore, the first law of thermodynamics combined with energy conservation tell us that 
the co-moving entropy density a is a constant. Note that the above results f[2T]) - fr23l) are 
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independent of the volume. They are derived based on the following assumptions: (1) The 
energy density and pressure are functions of temperature only; (2) The entropy is extensive 
so that we can define entropy density; (3) The volume and the temperature are independent. 

In the discussion of the phantom thermodynamics, it was assumed that the temperature 
of DE is proportional to the horizon temperature, i.e., T = 6T4 33j. For 6 = 1, the DE is in 
thermal equilibrium with the horizon. We adopt this temperature of the Universe T = bT& 
to discuss the DE thermodynamical property. Substituting T into Eq. (121]) . using Eqs. ffTUj) 
and (123]) . we get 

a = £. = g ~ k l a2 (24) 

a 3 2bG{H 2 + k/a 2 y/ 2 ' 1 } 

Considering the Friedmann equation (flOl) . we see that for quintessence with w > —1, s > 0, 
while for phantom with w < —1, s < 0. For the flat universe k = 0, the solution to Eq. ( |24j) 
is 

a(t) = [2baG(t-Q} 1/3 , (25) 
where is an integration constant, or 



a(t) 



3d(i-t„) 2krG 



nl/3 



(26) 



3d 

where d 7^ is an integration constant which is the asymptotical Hubble parameter. Only 
Eq. (|26l) gives the acceleration of the Universe driven by the DE. If a = 0, Eq. (|26|) becomes 
the cosmological constant solution. Substituting Eq. (T26|) into Friedmann equations (l8l)- (|T0l) . 
we get the DE equation of state 

AboG 



P = ~P 



_ e -3d(t-U) 



(27) 



3d 

Eq. (1271) is the DE equation of state when the DE temperature T = 6T4. The DE is in 
thermal equilibrium with the apparent horizon when 6 = 1. We see that when the DE is 
the cosmological constant, a = and p = —p. If the DE persists with an equation of state 
different from that described in Eq. (1271) . it concludes that DE is not in thermal equilibrium 
with the apparent horizon in general. 
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A. Dark energy with constant w 



For a fluid with constant equation of state p = wp, the solution to the energy conservation 
equation is p w — p w o(a / 'a) 3 ^ 1+w ^ . From Eq. (1221) we learn that 



, V(l+w) / \ -3 



S w \PwO / \ a 0/ 

where the subscript "0" denotes the present time value. The entropy inside the apparent 

horizon is g ^ 2 

^ = Us w f s A = f-) f— ) , (29) 
3 V a o/ \PtoJ 

where pt is the total energy density and the Friedmann equation (jSJ) was used to derive the 
last equation. During the radiation domination (RD), p t = p r ~ a~ 4 , Sr = SR Q~ 3 ^ 2 (a/a ) 3 
and 5*a = nr 2 A /G = S , J 4 ^7o 1 ( a / a o) 4 ) where fl r0 = p r o/Pto- During the matter domination 
(MD), p t = p m ~ cr 3 , so = S R0 £l^ /2 (a/a ) 3/2 and 5 A = S A o^mo( a / a o) 3 - During 
the w-fluid domination, p t = p w ~ a~( 1+,u ), then S# = Si?o^o^ 2 ( a / a o) 3<1+3w ^ 2 an d = 
SAoQ-Ka/ao) 3 ^^. During RD and MD eras, the entropy increases as the Universe expands, 
but when the DE dominates w < —1/3, the total entropy decreases as the Universe expands. 
It is easy to check that the entropy in each epoch of the evolution of the Universe is bounded 
by the apparent horizon entropy, once the entropy bound is satisfied at an earlier era. At 
very early epoch, the entropy bound is 

Sr - Sro q; 1/2 (^)<i. (so) 



OA J AO v a 

So a/a = T 70 /T 7 > S m n; o 1/2 /S Ao ~ 10~ 32 if we choose H = 72 Mpc" 1 km s" 1 , S m ~ 10 88 , 
Q r o ~ 10~ 4 and Sao ~ 10 122 , thus we have T 7 < 10 19 GeV. The relationship (128|) between s 
and p may be used to provide other useful information, for example, if there is an energy 
bound, then we can get an entropy bound, which was used to argue the minimum mass of 



a primordial black hole in 38]. 



Substitute Eq. ([28]) into Eq. (J2U), we get 

[l+w) / \ ~3w 

(31) 



T / n \ w /( 1 + w ) / n \ ~ 3w 

-L w I Pw \ I 



T w o \ Pwo J \ a o , 

where T w0 = (l + w)p w o/s w0 . For the radiation, w — 1/3, we get the familiar relation p ~ T 4 . 
For the DE, we see that it becomes hotter and hotter as the Universe expands. At present, 
the DE density p w0 is much greater than that of the radiation p r0 , the DE temperature 
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could be much higher than the current cosmic microwave background temperature unless w 
is very close to —1. 

In 28(, the authors find the spectral distribution E 1 /™ /[exp(E /T) ± 1]. Assuming that 



the DE is massless and follows either the Bose-Einstein or Fermi-Dirac statistics, we get 

W r J 1 '' dE , (32) 



r 2tt 2 J exp[E(P)/T] ± 1 6n 2 w J exp[E(P)/T\ ± 1 

where P is the momentum and g is the internal degrees of freedom. The number density is 
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n 



P 2 dP gT l ' w rv}l w - l du / T \ 1/w 

' no - . (33) 



2tt 2 J exp(£/T)±l 6n 2 w J e u ± 1 u V^o^ 

Therefore, we get the dispersion relationship E(P) ~ P 3w . For radiation, we recover the 
usual dispersion relationship E = P. For DE, we see that the greater the momentum, the 
smaller the energy This is another strange property of DE. It could be that we cannot apply 
the above reasoning to the DE, i.e., the DE may not follow the usual statistics or it is not 
massless. 



B. The generalized Chaplygin gas 



For the GCG model, p c = -Ajp a c (a > 0), Eqs. dHJ and ([22} give us that 



,a+l 



C(p« +1 -A), 



T = C-V(l+a) ( 1 



A 



a/a+l 



(34) 
(35) 



The above relationships were derived in 32j. In terms of the variables w c o and p c o, we can 
write A = — w c0 p" +1 . Eqs. ( 1341 and ( 1351) can be re- written as 



ScO 



Pc = Pc0 



a+1 



1 + W c0 



Pc_ 

PcO 



1+a 



+ W, 



1 1 + w C Q ( T 



WcO 



w c0 V T c 



(l+a)/a 



cO 



-l/(l+a) 



(36) 
(37) 



where T c0 = (1 + Wc0 )°/(i+«)c-V(«+i) = (i + Wc0 )p c0 / Sc0 . 
The solution to the energy conservation equation is 



Pc = PcO 



-w c0 + (1 + w 



cO ) 



a \ 3 ( 1 +«)' 



(38) 
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This equation tells us that the GCG energy density decreases as the Universe expands 
and p c — ► p c o(— w c o) l ^ 1+a ^ as a — ► oo, i.e., the GCG extrapolates between a dust and a 
cosmological constant. From Eq. (137)) . we see that the temperature of the GCG starts from 
the maximum temperature T m = T (l + w c o)~ a ^ 1+a ^ and decreases to zero as the Universe 
expands. The temperature behavior of the GCG is different from that of the DE with 
constant w. The temperature of the GCG decreases instead of increasing. This tells us that 
the increase or decrease behavior of the temperature of the Universe dominated by DE is 
model dependent, it is not the general property associated with the DE. 



IV. VOLUME AND TEMPERATURE INSIDE THE APPARENT HORIZON 

In cosmology, the physical volume is a function of the scale factor a(t). Since the scale 
factor a(t) relates to the temperature T of the Universe, the volume V should be a function 
of T. The first law ( 1191) becomes 

(Ts-p-p)dV = V(dp-Tds). (39) 

The formula (12TJ) - (|2"31) are special solution to Eq. (131?j) and they are not valid in general. As 
discussed in section HT1 w e should use the apparent horizon in the discussion of thermody- 
namics in cosmology 271 ] . For the radiation entropy, the previous results still hold. During 
RD, the entropy contributed by relativistic particles is Sr/Srq = Q r ^ 2 (a/a ) 3 . During MD, 
S R /S m = fi m o /2 (a/a ) 3/2 . Substituting V = Anr 3 A /3 into Eq. ([39]) and using Eq. ©, we get 

TdS = — - ( — - p7 5/2 [3(p + p)dp t - 2p t dp] 

, ( 3/2 ( 4 °) 
Pt {pt + 3p t )dp, 



3 V 3 

where the second equality was obtained by using the energy conservation law. During DE 
domination, TdS > if p < —p and TdS < if p > —p. 



A. Dark energy in thermal equilibrium 

If we think that DE is in thermal equilibrium with the Hawking radiation of the apparent 
horizon, i.e., T = Ta = I/Ikta, then we get 

dS=—^(p t + 3 Pt )p; 3 dp. (41) 
12 



For dark energy p w = wp w , the solution to Eq. (l4"Ij) during RD is 



3 ( X +W ) o n-2n ( a 



5— 3w 

S w = -z — - — -Sao^to^wo ( — ) • (42) 
5 — 3w \a J 

This is the DE entropy during RD period. During MD, the DE entropy can be obtained 
from Eq. (j4"Tl) in the form 

1 + W ( (X \ W ^ 

S w = S w \ + — -SAO^mO^wO I — I , (43) 

2(1 —w) \a J 

where S w \ is an integration constant. We see that the entropy of the DE increases faster 
than those of the radiation and the apparent horizon during RD and MD. To require that 
the current DE entropy is negligible compared to the radiation entropy, we get 

1i / \ 3/2-3™ a l/2 

l±w_ ( <<: „ 10 -34 ; (44) 



2(1 -w) \a J ~ Sao^ w o 

where a t is taken to be the end of the MD which is roughly the transition from deceleration 
to acceleration. The requirement of neglecting the DE entropy puts the constraint on the 
DE equation of state to be very close to —1. This is consistent with recent observations. 
If the DE equation of state is different from —1, then the DE entropy dominates the en- 
tropy components when the DE is in the thermal equilibrium with the apparent horizon. 
Comparing the DE entropy with the apparent horizon entropy 



S u 
5" 



a 



l + w Q w0 ( a t 



-3w 



~ < 1, (45) 



a~a t 2(1 — w)£l m0 \a 

we find that the holographic bound does not put strict constraint on w and it can be satisfied 
for general DE equation of state. 

During DE domination, the radiation entropy contributed by relativistic particles be- 
comes Sr/Srq = ^o^ 2 ( a / a o) 3 ^ 1+3M '^ 2 ) the apparent horizon entropy becomes Sa/Sao = 
^wl( a / a o) 3{1+w) and the solution to Eq. (jUJ) gives the DE entropy 

Sw = S wl + — (1 + 3w)p- 1 = S wt + -I^SaoQ-1 (-J , (46) 

where S W i is an integration constant. For DE with w < —1/3, the radiation entropy decreases 
while both the apparent horizon entropy and the absolute value of the DE entropy increase 
at the same rate. Eventually the radiation entropy will vanish and the DE entropy will 
become negative whose statistical meaning becomes hard to understand. The total entropy 

13 



of the space-time is S t = S w + Sa = S wi + 3(1 + w)SAo^Zo( a / a o) 3 ^ 1+w ^ /2- Therefore the 
total entropy is positive if w > — 1. To see the above results more clearly, we solve Eq. (14ip 
numerically for w = —0.9 and the result is shown in Fig. [IJ It is evident that S w increases 
faster than Sr and S w > Sr during RD and MD, and S w is negative after MD. The entropy 
bound is always satisfied. 

o 
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FIG. 1: The evolutions of S w , Sr and Sa- The dotted line is for \S w \/Sao (note that S w is negative 
after the MD), the dash-dot line is for Sr/Sao, the dashed line is for Sa/Sao, and the solid line is 
for the total entropy (Sa + Sr + S w )/Sao- 

For the GCG, Eq. (JHSJ) tells us that p c = p c0 (l + w c0 ) 1/(1+o) ( a o/a) 3 during RD and the 
solution to Eq. (f4T!) is 

S c = l&uM^l + ^ c o) 1/{1+Q) (£\ ° ■ (47) 
The entropy of the GCG increases faster than that of the radiation. 
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During the epoch of the GCG domination, the solution to Eq. (Pill) is 



S c — S c i + -Sao ( — — 
2 \pc 



Pc_ \ a+l 3w c0 



p c0 J a + 2 

S ci + IS A0 {1 + w c0 )- 1 /(i+«)(o/a ) 3 , a < a , 
^)(-w c o)- 1/(1+q) 5ao, a»ao, 



(48) 



where S ci is an integration constant. It is easy to see that S c0 = SAo(a+2+3w c0 )/2(a+2) and 
the holographic entropy bound is satisfied. If a > 1, then the entropy is always non-negative 
and it increases first and then decreases. The apparent horizon entropy is 



Sa — S 



AO 



-w c0 + (1 + w c0 ) 



«0 



-3(l+a) 



(1 + Wc0 )- 1 /( 1+Q ^ AO (a/ao) 3 , a < a , 
(-^ c o)- 1/(1+a) ^o, a»ao- 



(49) 



The entropy of the GCG is always smaller than that of the apparent horizon and the holo- 
graphic bound is always respected. The radiation entropy is 

-3/2(l+a) 

/ (I \ 

-w c0 + (1 + w c0 ) 



s R = s 



BO 



ao 



(1 + w c0 )-W^S r0 (a/a )^, a « a , 



(50) 



-w c0 )- 3 /^S r0 (a/a y 



The radiation entropy increases first and then decreases to zero. Since the entropy of the 
GCG increases faster than that of the radiation, the GCG entropy is greater than the 
radiation entropy at the present because S c0 > S r0 and the holographic entropy bound is 
always respected. These results are evident from Fig. [2] with w c0 = —0.88 and a = 1.57 



34|. So the GCG entropy is the dominate entropy if it is in thermal equilibrium with the 



apparent horizon. 



B. Dark energy with power law entropy 



Recall that the radiation temperature T ~ a 1 , the radiation entropy S ~ T 3 during 
RD and S ~ T~ 3 / 2 during MD. The relationship between the entropy and the temperature 
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FIG. 2: The evolutions of S c , Sr and Sa with w c q = —0.88 and a = 1.57. The dotted line is for 
Sc/Sao, the dash-dot line is for Sr/Sao, and the solid line is for Sa/Sao- 



is a power law form. More generally, we assume that S = (3T 1 ', where the constants j3 and 
7 depend on the background, then Eq. ( HOI) becomes 



-3/2 



For the DE with p = wp, the solution to Eq. ( IST1) during RD is 



w 1 



4tt(1 +w) (MG 



-3/2 x x 3(1— to) 

-3/2 / a 

PrO PwO I 

V«0, 



(51) 



(52) 



7^1 + 1"™ 3(1 — w) V 3 

For radiation, j w i = —3 and T w ~ a -1 . If j w i < —1 for —1/3 > w > —1, then the DE 
temperature decreases and the entropy increases as the Universe expands. Since Sr ~ a 3 and 
~ a 3(i-w)7 lt ,i/(i+7 m i) > fl 3^ j-j-g en tropy increases faster than the radiation entropy. 

To make sure that the DE entropy is negligible compared to the radiation entropy, we 
require < j w i < —1/w so that j w i/(l + j w i) > and 3(1 — w)j w i/(l + 7^1) < 3 
for —1/3 > w > —1. Therefore the DE temperature and entropy both increase as the 
Universe expands. For w < —1, we need — 1 < 7^1 < to keep the temperature positive. 
The phantom temperature increases and the entropy decreases as the Universe expands. 
However, both the temperature and the entropy of the phantom are positive. 
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During MD, the solution to Eq. (|5Tj) is 



7^2 + 1 3(1 — 2u>) \ 3 J \a ) 

where C\ is an integration constant. If 7^2 < — 1 for — 1/3 > w > —1, then the DE entropy 
increases at a faster rate than that of the radiation as the Universe expands. So we require 
< 7^2 < -l/2w to keep 7^2/(1 + 7^2) > and 3(1 - 2w)j w2 /2(l +j w2 ) < 3/2. Again both 
the temperature and entropy of the DE with w > — 1 increase as the Universe expands. For 
w < —1, we require —1 < 7^,2 < 0. The phantom temperature increase and the entropy 
decreases as the Universe expands. 

During the DE domination, the solution to Eq. ( 15TT) is 



7^,3 + 1 3 V 3 / V a o/ 

where C2 is an integration constant. Now we require that —1 < 7^3 < for both the 
quintessence and phantom to get a positive temperature. For w > —1, the temperature 
increases and the entropy decreases as the Universe expands. For w < — 1, the temperature 
decreases and the entropy increases as the Universe expands. At late times, the entropy is 



u>3 

a 
a 



I /o / <\-3/2 . , / \ 3(l+™)/2l 7»s/(7»3+l) 

4vr / 87tG\ 7^3 + 1 / a x 

l + 3w)- 



3/?to3 \ 3 J 7^3 \ao 

3(l+w) 7 „ 3 /2(l+7»3) 



(55) 



Even for the phantom, both the temperature and the entropy are positive. For w > — 1, 
during the DE domination, the radiation entropy decreases as a 3 ( 1 + 3w )/ 2 ) the DE entropy 
decreases as a 3< - 1+1u ^ 7//2 ^ 1+7 ^ and the horizon entropy increases as a 3( - 1+w \ so the entropy 
bound will be respected. As we discussed earlier, the entropies of the radiation and DE with 
w > — 1 increases slower than that of the apparent horizon, so the entropy bound is always 
satisfied. While for the phantom, w < — 1, both the horizon entropy and the radiation 
entropy decrease and the phantom entropy increases, the entropy bound will be violated at 
late times. 

For the GCG, the solution to Eq. (15T]) during RD is 

^' +1 -f(¥)" 3/2 ( 1+ -) w+ " , - / -fe) 3 - («) 
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So the temperature scales as T c ~ a 3 ^ 1+7c1 ^ and the entropy scales as S c ~ a 3 7ci/(i+7ci)_ 
If 7ci < —1, then the entropy increases faster than that of the radiation as the Universe 
expands. If J& > 0, then both the temperature and the entropy increases and the entropy 
increases slower than that of the radiation as the Universe expands. During the GCG 
domination, the solution to Eq. (1511) is 



^« _o. + * (?f)- 3/ > {ff + a + ^ ( 

(57) 

where C3 is an integration constant. At early times, a C 00, the temperature scales as 
T c ~ a 3 / 2 ( 1+7c2 ) and the entropy scales as S c ~ a 37c2 / 2 ( 1+7c2 ). At late times, a > a , both 
the temperature and the entropy approach to a constant. To require that the DE entropy is 
negligible compared to that of radiation at least up to the present, we find that 7 C 2 > 0. So 
both the DE temperature and entropy increase first and then decrease and the DE entropy 
increases slower than that of the radiation. 



'a \ 3 ( 1+Q ) 
a 



V. THE GENERALIZED SECOND LAW 

Now we proceed to discuss the GSL for the Universe with DE. The geometric entropy 
associated with the apparent horizon is Sa = nr\/G. Using Eq. ( fTTl) . we get 

S A = 3S A H^^. (58) 
Pt 

As long as p t + Pt > 0, Sa > , while for the phantom domination, Sa < 0. 



A. Constant co-moving entropy density 

The radiation entropy contributed by the relativistic particles inside the apparent horizon 
is S = 4iTsr\/3, its time derivative reads 

S = SSH^P^. (59) 
2pt V ; 

In the radiation or matter dominated eras, S increases as the Universe expands if S > 0, so 
the GSL is satisfied. The violation of the GSL is only possible when DE dominates. During 
the DE domination, the variation of the total entropy in the Universe is expressed as 

S + S A = 3SH*±^ + 3S A H^±^. (60) 
2pt pt 
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During RD and MD eras, the radiation entropy increases slower than that of the apparent 
horizon, so the entropy bound S < Sa is always satisfied. Considering the equation of state 
in the RD and MD epoches, Eq. (IBDj) is always positive so that the GSL is protected during 
RD and MD eras. When the DE starts to dominate, S decreases as the Universe expands 
while Sa keeps on increasing for w > —1, so in the far future S <C Sa for w > — 1. Since 
(pt + 3pt)/2pt has the same order of magnitude as (p t + Pt)/pt, the GSL can be preserved 
for w > —1. To be more explicit, we take the DE model with constant w as an example. In 
this case, Eq. (I60I) becomes 



S + S A = 3H 



l + 3w Q /a\ 3(1+3w)/2 n ^ fa^ 3il+w) ' 
~^^S — + (1 + w)S A0 - 

2 \a / \a 



(61) 



For w > — 1, |(1 + 3w)/2\S < (1 + w)Sao an d the first term in the right hand side decreases 
and the second term increases, so the GSL is always satisfied. 

For the phantom with w < —1, from the discussion in Sec. IIII A\ we learnt that the 
temperature T < while the entropy S > 0. The GSL is violated since Sa < and 
5 < during the phantom domination. This problem can be attributed to the negative 
temperature deduced in the formalism where the volume and the temperature are assumed to 



be independent. In 



33| , by identifying the temperature to the apparent horizon temperature, 



T > and S < 0, it was shown that the GSL may be satisfied. 



B. Volume as a function of Temperature 

Now we study the GSL by considering the physical volume and the temperature are 
dependent and taking the temperature of the Universe to be equal to the apparent horizon 
temperature. From Eq. (1401) . we see that S > (< 0) if (p + p)/T > (< 0) during RD 
and MD, and S < (> 0) if (p + p)/T > (< 0) during the DE domination. During the 
phantom domination, Sa < 0, so we must require T > to protect the GSL. 

If T = T A , Eq. flUD tells us that 

S= Y^(pt + 3 Pt )(p + p)pT 3 H. (62) 
Combining Eqs. fl58l) and fl62l . we get 
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St + S A = j^Hpfipt + Pt ) 2 > (63) 
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for the total entropy of the Universe and the apparent horizon entropy. Therefore, the GSL 
is always satisfied even for the phantom. 

For the case of S — (3T y , we find that the entropies all increase during RD and MD. 
Only when the DE dominates, the entropies of the radiation and DE with w > — 1 decrease. 
So we discuss the GSL during the DE domination. All the entropies have the power law 
dependence on the scale factor a, so S ~ SH with an order one proportionality factor. Since 
Sr <C Sa and S w <C Sa, so the GSL is always respected. 

VI. CONCLUSIONS AND DISCUSSION 

In this paper we have studied the thermodynamical properties of the Universe with DE. 
Adopting the usual assumption in deriving the constant co-moving entropy density that 
the physical volume and the temperature are independent, we studied the thermodynamics 
of the Universe in the radiation, matter and DE dominated periods. During RD and MD 
epoches, the entropy of the Universe contributed by the relativistic particles increase with 
the expansion, while in the DE dominated era, the entropy remains positive but decreases 
as the Universe expands. In all these processes, the entropy within the Universe is bounded 
by the geometric entropy associated with the apparent horizon. The holographic principle 
is respected. The requirement of the holography at early epoch gives reasonable energy 
scale. We have found some strange behaviors of the DE based on the assumption that 
the physical volume and the temperature are independent. For the DE with constant w, 
we have observed the strange dispersion relation E ~ P 3w . For the phantom, we saw 
that either the temperature or the entropy could be negative. For the DE with constant 
equation of state w > — 1, the temperature increases as the Universe expands. We have also 
extended our discussion to the GCG and found that its temperature decreases as the Universe 
expands, which is different from the DE with constant w. Thus the strange behavior on the 
temperature as observed is not the general property of the DE, it looks model dependent. 

We have also considered the realistic situation that the physical volume and the temper- 
ature of the Universe are related. We have concentrated on the volume within the apparent 
horizon. The apparent horizon is important for the study of cosmology, since on the appar- 
ent horizon there is the well known correspondence between the first law of thermodynamics 
and Einstein equation. On the other hand it has been found that the apparent horizon is 
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a good boundary for keeping thermodynamical laws [27]]. Considering that the Universe is 
in thermal equilibrium, T = Ta, we have studied the entropy of the DE in the radiation, 
matter and DE dominated eras. We found the DE entropy increases during the RD and 
MD eras while decreases during the DE dominated era. In all epoches, the DE entropy is 
bounded by the holographic entropy on the apparent horizon. Requiring that at the present 
epoch, the DE entropy is negligible compared to the radiation entropy, we can limit the DE 
equation of state to be very close to —1, which is consistent with observations. We found 
that both temperature and entropy can be positive for the DE including the phantom. 

We have also investigated the GSL in the Universe with DE. By assuming the physical 
volume independent of the temperature, we have found that the GSL is protected for the DE 
with w > — 1. Considering the realistic case that the physical volume and the temperature 
are related, and the DE inside the Universe is in thermal equilibrium with the boundary by 
the apparent horizon, the GSL is proved to be always satisfied. 

The results tell us that in studying the thermodynamics of the Universe with DE, it is 
more appropriate to consider the case that the physical volume and the temperature are 
related. This could give us more reasonable results on the thermodynamical quantities in 
the Universe with DE. 
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